Orbital dynamics in ferromagnetic transition metal oxides 



00 
On 
OS 

o 

Q 

00 



3 



c 

o 
o 



> ■ 
m ■ 

<n : 

(N ■ 

oo : 

Os . 

-i— > ■ 

a ■ 

a 

T3 ■ 

c : 
o . 
o ■ 



X 



Jeroen van den Brink, Peter Horsch, Frank Mack, and Andrzej M. OlesB 
Max-Planck-Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Federal Republic of Germany 

(September 22, 1998) 

We consider a model ol strongly correlated e 9 electrons interacting by superexchange orbital 
interactions in the ferromagnetic phase of LaMn03. It is found that the classical orbital order 
with alternating occupied e 9 orbitals has a full rotational symmetry at orbital degeneracy, and the 
excitation spectrum derived using the linear spin-wave theory is gapless. The quantum (fluctuation) 
corrections to the order parameter and to the ground state energy restore the cubic symmetry of the 
model. By applying a uniaxial pressure orbital degeneracy is lifted in a tetragonal field and one finds 
an orbital-flop phase with a gap in the excitation spectrum. In two dimensions the classical order 
is more robust near the orbital degeneracy point and quantum effects are suppressed. The orbital 
excitation spectra obtained using finite temperature diagonalization of two-dimensional clusters 
consist of a quasiparticle accompanied by satellite structures. The orbital waves found within the 
linear spin-wave theory provide an excellent description of the dominant pole of these spectra. 

PACS numbers: 71.27.+a, 75.10.Jm, 75.30.Ds, 75.40.Gb. 



I. INTRODUCTION 

Recently much attention has been attracted to the 
properties of Lai-^A^MnOs, with A=Ca,Sr,Ba, and re- 
lated compounds, in which the colossal magnetaresis- 
tance and metal-insulator transition are observed.EI The 
parent compound, LaMn03, is insulating with layered- 
type antiferromagnetic (AF) ordered state (the A-type 
antiferromagnet), accompanied by the orbital order with 
alternation of fig orbitals occupied by a single electron 
of Mn 3+ ions.B Such orbitally ordered states pBiight be 
promoted by a cooperative Jahn- Teller effect Jj playing 
a decisive role in charge transportJj but could also re- 
sult from an electronic instability!! The latter possibil- 
ity is now under debate, as the band structure calcu- 
lations performed within the extension of local density 
approximation (LDA) to strongly correlated transition 
metal oxides, the so-called LDA+U approach, give in- 
deed charge-ordered ground states without assuming lat- 
tice distortions as their driving force. □ Charge ordering 
was also found within the Hartree-Fock calculations on 
lattice models by Mizokawa and Fujimori.Ll 

The dominating energy scale in late transition metal 
oxides is the local Coulomb interaction U between d 
electrons at transition metal ions which motivates the 
description of ground state and low-energy excitations 
in terms of effective models with spin and orbital de- 
grees of freedom. Such models, introduced both fop-jthe 
cuprates Jao and more recently for the manganites ,II2lil2l 
show that the spin and orbital degrees of freedom are 
interrelated which leads to an interesting problem, even 
without doping. Orbital interactions in the ground state 
lead to a particular orbital ordering that depends on the 
magnetic ordering, and vice versa. An interesting aspect 
of such system is that the elementary excitations have 
either a puceuspin, or pure orbital, or mixed spin-orbital 
char act er.Ej~li3 This gives large quantum fluctuation cor- 
rections to the order parameter in AF phases, and it was 
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Mi 



discussed recent 

of a spin-liquid£E9'[Li3 In a special case of ferromagnetic 
(FM) order (which represents a simplified spinless prob- 
lem) , only orbital excitations contribute to the properties 
of the ground state and the quantum effects are expected 
to be smaller. 

Unlike in the cuprate materials, hole doping in the 
manganites leads to FM states, insulating at low dop- 
ing, and metallic at higher doping. E_ZJ Although the FM 
phase is unstable for the undoped LaMnC>3, the build- 
ing blocks are FM planes. A proper understanding of 
pure orbital excitations has to start from the uniform 
FM phase, where the spin operators can be integrated 
out. Such a state could possibly be realized in LaMn03 
at large magnetic field and will serve as a reference state 
to consider hole propagation in doped manganites, simi- 
larly as the AF spin state provides a reference state for 
the hole propagation in the cuprates. Therefore, we con- 
sider in this paper an effective model with only orbital su- 
perexchange interactions which involves the e g states and 
results from the effective Hamiltonians derived in Refs. 
p|-p^ in the relevant high-spin state. Of course, it is iden- 
tical for the cuprates and for the manganites. We note, 
however, that orbital interactions may be also induced by 
lattice distortions, but these contributions are expected, 
to be less important than the orbital superexchange ,113 
and will be neglected. 

So far, little is known about the consequences of orbital 
excitations for long-range order and transport properties 
in doped systems, but it might be expected that such 
excitations either bind to a moving hole, ox_lead to hole 
scattering even when the spins are aligned.Ea The knowl- 
edge of orbital excitations is a prerequisite for a better 
understanding of the temperature dependence of the re- 
sistivity and optical conductivity in the FM phase of 
doped manganites. Although the double exchange mech- 
anism is quite successful in explaining why the spin order 
becomes FM with increasing doping of Lai_ ;E A. c Mn03, it 
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fails to reproduce the experimentally observed tempera- 
ture dependence of the resistivity] 2 !] Even more puzzling 
and contradicting naive expectations is the incoherent 
optical conductivity with a small Drude -peak observed 
in a FM metallic phase at higher doping.Eil This behav- 
ior demonstrates the importance of orbital dynamics in 
doped manganites which appears to be predominantly 
incoherent in the relevant parameter regime 

The paper is organized as follows. The orbital superex- 
change model is presented in Sec. II. We derive the clas- 
sical phases with alternating orbitals on two sublattices 
in three and two dimensions, and analyze their depen- 
dence on the crystal field splitting. Orbital excitations 
are derived in Sec. Ill using a pseudospin Hamiltonian 
and an extension of the linear spin- wave theory (LSWT) 
to the present situation with no conservation of orbital 
quantum number. We present numerical results for the 
dispersion of orbital waves and determine quantum fluc- 
tuation corrections to the ground-state energy and to the 
order parameter in Sec. IV. The derived spectra are com- 
pared with the results of exact diagonalization in Sec. V. 
A short summary and conclusions are given in Sec. VI. 



II. ORBITAL HAMILTONIAN AND CLASSICAL 
STATES 

A. Orbital superexchange interactions 

We consider a three-dimensional (3D) Mott insulator 
with one e g electron per site. The effective Hamiltonian 
which descpbes e g electrons in a cubic crystal at strong- 
coupling Jj~ll3 

h = Hj + h z , (i) 

consists of the superexchange part Hj, and the orbital 
splitting term due to crystal-field term H z . Here we will 
consider a special case of the FM e g band in which the 
spin dynamics is integrated out and one is left with the 
effective interactions between the electrons (holes) in dif- 
ferent orbital states. In this case the only superexchange 
channel which contributes is the effective interaction via 
the high-spin state, and thus the interaction occurs only 
between the pairs of ions with singly occupied orthogo- 
nal e g orbitals at two nearest-neighbor sites (two orthog- 
onal e g orbitals are then singly occupied in the interme- 
diate excited states). An example of such a model is the 
superexchange interaction in the FM state of LaMnC>3 
which originates from djdj dfd^ excitations (with 

d t = 4g e g^ d ! = 4g, and d t = 4 g e 2 g ) into a high-spin 
d| state | 6 Ai), as shown schematically in Fig. |]. It eives 
the effective Hamiltonian with orbital interactionsJ3'El 



where t is the hopping element between the directional 
3z 2 — r 2 orbitals along the considered c-axis, and £( 6 ^4i) is 
the excitation energy. The orbital degrees of freedom are 
described by the projection operators 'P/fj) which select a 
pair of orbitals |£) and being parallel and orthogonal 
to the directions of the considered bond (ij) in a cubic 
lattice. 

The Hamiltonian ^ has cubic symmetry and may be 
written using any reference basis in the e g subspace. For 
the conventional choice of 3z 2 — r 2 ~ \z) &ndx 2 — y 2 ~ \x) 
orbitals, the above projection operators are represented 
by the orbital operators r", with a = a,b,c for three 
cubic axes, 

?%) = i\- r?)G + r?) + (§ + rf )(§ - rf). (3) 

It is convenient to replace the orbital operators rf by the 
pseudospin operators Tf with ji = x,z, 

rf b) TN/3T?), if = T i ', (4) 

The latter operators may be represented by the Pauli ma- 
trices in the same way as spin operators, Tf = \(rf 1 and 
obey the same commutation relations. Therefore, the zth 
component of pseudospin is given by Tf — \{ni X — rii Z ), 
and we identify the orbital states as up- and down- 
pseudospin, \x) = | T)i an d \z) = | j), respectively. We 
take the prefactor J = t 2 /e( 6 Ai) in Eq. (||) as the energy 
unit for the superexchange interaction. Thus, one finds 
a pseudospin Hamiltonian, 

Hj = | J [ T " T I + ?> T i T i T VE(T?Tf + T?Tf) 
«>ll 

+ 2J y, T t T h ( 5 ) 

(ij)-L 

where the prefactor of the mixed term cx \[2> is neaa-, 
tive in the a-direction and positive in the 6-direction.Efl 
We choose a convention that the bonds labeled as 
((ij)A.) connect nearest-neighbor sites within (a, b) planes 
(along the c-axis). The virtual excitations which lead to 
the interactions described by Eq. (||) are shown in Fig. 
|l|. Here we neglected a trivial constant term which gives 
the energy of —J/2 per bond, i.e., — 3J/2 per site in a 
3D system. We emphasize that the SU(2) symmetry is 
explicitly broken in Hj, and the interaction depends only 
on two pseudospin operators, Tf and Tf . 

The crystal-field term removes the degeneracy of \x) 
and \z) orbitals, 

H z = -E z J2Tf, (6) 

i 

and is connected with the uniaxial pressure. This term is 
induced by static distortions in a tetragonal field and is of 
particular importance in a two-dimensional (2D) system. 
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B. Classical ground states 

Before analyzing the excitation spectra of the Hamil- 
tonian (||) , one has to determine first the classical ground 
state of the system that will serve as a reference state to 
calculate the Gaussian fluctuations. Neglecting the ir- 
relevant phase factors, the classical configurations which 
minimize the interaction terms (^|) are characterized by 
the two-sublattice pseudospin order, with two angles de- 
scribing orientations of pseudospins, one at each sublat- 
tice. As usually, the classical ground state is obtained by 
minimizing the energy with respect to these two rotation 
angles, i.e., by choosing the optimal orbitals. 

Let us consider first the term Hj at orbital degeneracy 
E z = 0. The superexchange interaction (|J) induces the 
alternation of orthogonal orbitals in the ground state in 
all three directions, and is equivalent to the two-sublattice 
AF (G-AF) order in the pseudospin space. This config- 
uration gives the lowest energy on the mean-field level, 
as the virtual transitions represented in Fig. 1 give the 
largest contribution, if the hopping involves one occupied 
and one unoccupied orbital of the same type (e.g., either 
directional or planar). A priori, the optimal choice of 
the occupied orbitals could be unique on the classical 
level due to the cubic symmetry of superexchange inter- 
action (Q). Therefore, we perform a uniform rotation of 
{\z}, \x}} orbitals by an angle at each site, 
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to generate the new orthogonal orbitals, \ip) and \iv), 
and investigate the energy as a function of 0. The AF 
order imposed by the Hamiltonian Hj in a cubic system 
implies the alternation of these orthogonal orbitals in the 
rotated basis, \ip) and \ju), at two sublattices, i.e., i G A 
and j € B. A few examples of various possible choices 
of the orthogonal orbitals {\ip), \iv)} are given in Table 
I for representative values of 0. 

The rotation (]?]) leads to the following transformation 
of the pseudospin operators, 



Tf — > Tf cos 28 — T* sin 28, 
Tf -► Tf sin 20 + Tf cos 20, 



(8) 



and the interaction Hamiltonian Hj is then transformed 
into, 



H e ± = jJ2 i i 1 ~ cos40 ) T i T j + i 1 + cos46) ) T i T j 
(ij)-i- 



-sm40(TfT/ + lfT?) 



(11) 



n° = m + m 



(9) 



As in Eq. (Js[) , the bonds in the first sum (|ic| ) are par- 
allel to either a or 6-axis, while the bonds in the second 
sum (pd| ) are parallel to the c-axis. The classical energy 
is minimized if (TfT?) — — |, i.e., if the orbitals order 
'antiferromagnetically'. 

The Hamiltonian given now by Eqs. (fuj and @ has 
the symmetry of the cubic lattice, but surprisingly one 
finds the full rotational symmetry of the present interact- 
ing problem on the classical level at orbital degeneracy. 
In other words, classically the lowest energy is always 
Emf — —3 J/4 per site, independent of the rotation an- 
gle, as long as alternating orbitals on neighboring sites 
are occupied. This follows from the particular structure 
of the rotated Hamiltonian @ which has identical factors 
of three in front of the diagonal cx T?T? and off-diagonal 
oc T?Tj contributions, when these are summed over all 
the bonds which originate at each site, and these coeffi- 
cients are independent of the rotation angle 8. In contrast 
to the Heisenberg antiferromagnet (HAF), however, this 
symmetry concerns only the summed contributions, and 
not the interactions along the individual bonds. 

A finite orbital field E z ^ breaks the rotational sym- 
metry on the classical level. It acts along the c-axis, and 
it is therefore easy to show that the ground state in the 
limit of E z — ► is realized by the alternating occupied 
orbitals being symmetric/antisymmetric linear combina- 
tions of \z) and \x) orbitals,E3 i.e., the occupied states 
correspond to the rotated orbitals (|7]) \ip) and \iv) on 
the two sublattices with an angle 8 = 7r/4, shown in Fig. 
§ (see also Table I). In particular, this state is different 
from the alternating directional orbitals, 3x 2 — r 2 and 
3y 2 — r 2 , which might have been naively expected. It fol- 
lows in the limit of degenerate orbitals from the 'orbital- 
flop' phase, in analogy to a spin-flop phase for the HAF 
at finite magnetic field. With increasing (decreasing) E z 
the orbitals tilt out of the state shown in Fig. 0, and 
approach \x) (\z)) orbitals, respectively, which may be 
interpreted as an increasing FM component of the mag- 
netization in the pseudospin model. 

It is convenient to describe this tilting of pseudospins 
due to the crystal-field oc E z by making two different 
transformations (0) at both sublattices, rotating the or- 
bitals by an angle 8 — j — <j> on sublattice A, and by an 
angle = ? + <j> on sublattice B, so that the relative an- 
gle between the occupied orbitals is ? — 20 and decreases 
with increasing 4>, i.e., with increasing E z , 



Hj = - K 2 + 00846 T V3 sin40)Tf T/ 

(«>ll 

+ (2 - cos 40 ± v / 3sin40)Tf T? 
-(sin 40 ± V3 cos 48) (TfT* + T?Tf) 



(10) 



MW cos(f-. 

-sin(f - 



\jy) 



cos(f + ( 
-sin(| + , 



sin(f 
cos(f 



sm(f 
cos(f 



\IZ) 
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(12) 



(13) 
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As before, the orbitals and \jv) are occupied at two 
sublattices, i £ A and j S B, respectively, and the orbital 
order is AF in the classical state, with the transformed 
operators (2?) = -1/2 and (T?) = +1/2, respectively. 
As a result one finds that the angles for the occupied 
orbitals are indeed opposite on the two sublattices. 

The operators 7? and Tf may be now transformed as 
in Eqs. (||) using the actual rotations by 9 — j ± 0, as 
given in Eqs. ( |l2| ) and (|l3|), and the transformed Hamil- 
tonian takes the form, 

— Ht + Hf_ + Hf, (14) 



H \\ = 2 E [(2 cos 40 - l)TfTj + (2cos40 + l)2?T; 

(«>ll 

+2sm40(TfT?-2?T/) 

iV3(Tf T/ + 2?T?) ] , (15) 

H ± = 3 E t ( cos4< ^ + i)^/ + ( cos4 <^ - i)^ 1 ; 2 

- sin 40(2? T? - 2?T?) ] , (16) 
Hf = E z J2( X i sin20Tf - cos20Tf), (17) 

i 

where \ = —1 for i G A and A, = 1 for i E B. The 
energy of the classical ground state is given by, 

££r7 = -^Jcos40-^sin20, (18) 

and is minimized by, 

sin 20= gj. (19) 

At orbital degeneracy (E z = 0) this result is equivalent 
to cos 2$ = [6 = 7r/4) in Eq. (Q), as obtained also in 
Ref. O, and one recovers from Eq. ( |l8| ) the energy of 
—3 J/4 as a particular realization of the degenerate clas- 
sical phases with alternating orthogonal orbitals. The 
above result (|l9|) is valid for \E Z \ < 6 J; otherwise one of 
the initial orbitals (either \x) or \z)) is occupied at each 
site, and the state is fully polarized (sin 20 = ±1). The 
case of = corresponds to the occupied orthogonal 
orbitals at orbital degeneracy. In contrast, the value of 
= 7r/4 corresponds to performing no rotation on the A 
sublattice, while the orbitals are interchanged on the B 
sublattice [see Eqs. ( |l2|) and (|l3|)]. The latter situation 
describes a state obtained in strong orbital field, with 
only one type of orbitals occupied, a 'FM orbital state'. 
The orbital field changes therefore the orbital ordering 
from the AF to FM one in a continuous fashion. 



C. Two-dimensional orbital model 

As a special case we consider also a 2D orbital model 
with the interactions in the (a, b) plane. In this case 
the cubic symmetry is explicitly broken, and the classical 
state is of a spin-flop type. It corresponds to alternat- 
ingly occupied orbitals on the two sublattices, with the 
orbitals oriented in rthe plane and given by 9 = tt/4 at 
E z = 0, see Fig. |[Ea A finite value of E z tilts the or- 
bitals out of the planar |x) orbitals by an angle 0, and 
the Hamiltonian reduces to 

Ht D = H{ + Hf, (20) 

as there is no bond in the c-direction. The classical en- 
ergy is 

B 2 M D F = -ij(2cos40+ 1) - sin 20. (21) 

Therefore, one finds the same energy of —3 J/4 as in a 3D 
model at orbital degeneracy. This shows that the orbital 
superexchange interactions are geometrically frustrated, 
and the bonds in the third dimension cannot lower the 
energy, but only allow for restoring the rotational sym- 
metry on the classical level and rotating the orthogonal 
orbitals in an arbitrary way. The energy ( |2l| ) is mini- 
mized by, 

sin20=|i, (22) 

if \E Z \ < 4 J; otherwise sin 20 = ±1. Interestingly, the 
value of the field at which the orbitals are fully polar- 
ized is reduced by one third from the value obtained in 
three dimensions (|l9|). This shows that although the 
orbital exchange energy can be gained in a 3D model 
only on the bonds along two directions in the alternating 
(orbital-flop) phase either at or close to E z = 0, one has 
to counteract the superexchange on the bonds in all three 
directions when the field is applied. 

We note that the present orbital model becomes purely 
classical in one dimension, if the lattice effects are ne- 
glected, as we have assumed in the present model. This 
follows from the derivation which gives the superex- 
change interactions of Ising type between the \z) orbitals 
along the c-axis, and one can always choose this axis 
along the considered chain. Thus, one finds the Hamil- 
tonian 

H 1D = 2JY, T * T *+1- E *Y, T *- ( 23 ) 

i i 

At orbital degeneracy (E z = 0) the occupied orbitals al- 
ternate along the one-dimensional (ID) chain between 
two orthogonal orbitals, and one may for example choose 
the occupied states as \z) and \x) on even and odd sites, 
respectively. The energy of this ground state per site is 
—J/2, which allows to conclude that energy gains due to 
orbital ordering per site are significantly reduced com- 
pared to the 2D and 3D case. 



4 



III. ORBITAL EXCITATIONS IN LINEAR 
SPIN- WAVE THEORY 

The superexchange in the orbital subspace is AF and 
one may map the orbital terms in the Hamiltonian (|^) 
onto a spin problem in order to treat the elementary 
excitations within the LSWT. The equations of motion 
can then be linearized by a standard technique.B Let 
us first discuss the results of LSWT for the spin-flop 
phase induced by an orbital-field, starting from the ro- 
tated Hamiltonian ([hi]), rlfere we choose the Holstein- 
Primakoff transformationcS for localized pseudospin op- 
erators (T =1/2), 



TP 



TP 



(1 



a\ai) 1/2 ai, 



rj^Z — T — 

T i =a\ a i- 2' 
for i £ A sublattice and 

T/ = (l-%) 1 /%, T~ 



(24) 



(25) 



for j £ B sublattice. After replacing the square roots 
by the leading lowest order terms, and inserting the ex- 
pansion into Eq. (fi4|), one diagonalizes the linearized 
Hamiltonian by two consecutive transformations. Note 
that, contrary to the HAF, also cubic terms in boson op- 
erators are present in the expansion due to the anomalous 
interactions ~ TfTf. 

The linearized Hamiltonian simplifies by introducing 
the Fourier transformed boson operators 5k and &k given 
by 



6w = 



jeB 



and by transforming to new boson operators {ak, 6k}, 

a k = (a k ~ 5k)/V2, &k = (a k + 6k)/V2. (26) 
One finds the effective orbital Hamiltonian of the form, 

#lsw = J ^ -4ka k a k + ^• B k(a k a -k + a k a -k) 



J 



k 

£ 

k 



2 



M-k) 



(27) 



where the coefficients A-^ and B^ depend on angle 0, 

(28) 

:os40 + l) 7± (k)] , (29) 



A k = 3 - B k , 

' cos 40 — l)7||(k) 



B 



[(2. 



and the k-dependence is given by 711 (k) = ^(cosfc^ + 
cosfcy), and by 7_i_(k) = cosfc z , respectively. After a Bo- 
goliubov transformation, 



a k = u k Q'k + «kaL k ) 
with the parameters 



(30) 



Uk 




- sgn( W^H' (31) 



with 77k = y A^ — B^ , and an equivalent transformation 
for the 6k bosons, the Hamiltonian ( p7| ) is diagonalized, 
and takes the following form, 



LSW 



k 



[w^(^)4«k+Wk(0)/3i/i]. (32) 



The orbital-wave dispersion is given by 

£^(0) = 3 J {1 ± i [(2 cos 40 - 1) 7 || (k) 
+ (cos40 + l) 7± (k)]} 1/2 . 



(33) 



The orbital excitation spectrum consists of two branches 
like, for instance, in an anisotropic Heisenberg model.EZI 
The dependence on the field E z is implicitly contained 
in the above relations via the angle 0, as determined by 
Eqs. © and @ for the 3D and 2D model, respec- 
tively. The dispersion changes between the case of = 
(cos 40 = 1) which stands for the spin- flop phase at or- 
bital degeneracy, and the case of = ir/4 (cos 40 = — 1) 
which corresponds to the uniform phase with either \x) 
{E z > 6J) or \z) (E z < — 6J) orbitals occupied. The or- 
bital field therefore changes the orbital-wave dispersion 
from AF to FM state. The orbital-wave dispersion for a 
2D system, can easily be found from Eq. ( |33| ) by setting 
7.1 (k) = 0. 

In the case of a vanishing orbital field the classical 
ground state is degenerate with respect to a rotation of 
the orbital around an arbitrary angle 9. In LSWT, nev- 
ertheless, the orbital- wave dispersion depends on 8, i.e., 
the orbital-wave velocity is highly anisotropic, and one 
finds, 



fi±(0) = 3j{l± i (V3cos26»-sin26») 2 7 x (k) 

(\/3cos20 + sin20) 2 7 y (k) + 4sin 2 20 7± (k) j 



+ 



(34) 



with 7x)9 (k) = cosfc^j, respectively. It is straightforward 
to verify that the excitations are gapless, independent of 
the rotation angle. If 9 = 0, the orbital-waves given by 
Eq. ( p3| ) take a particularly simple form, and is identical 
to the dispersion found in the spin- flop phase at = 7r/4, 



n±(0 = 0) = 3 JJl ± 7 || (k) = c^(0 = tt/4). 



(35) 



In both cases the system is quasi-2D and the dispersion 
originates only from the planar components k x and k y . 
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IV. NUMERICAL RESULTS 

The excitation spectra given by Eq. (|33| ) for different 
crystal-field splittings are shown in Fig. Hfor the 3D sjss-i 
tem along different directions of the fee Brillouin zoneEa 
appropriate for the alternating orbital order. Most in- 
terestingly, a gapless orbital-wave excitation is found for 
the 3D system at orbital degeneracy. Obviously this is 
due to the fact that the classical ground state energy is 
independent of the rotation angle 9 at E z — 0, as we 
have shown in Sec. H At first glance, however, one does 
not expect such a gapless mode, as the Hamiltonian (||) 
does not obey a continuous SU (2) symmetry. The cubic 
symmetry of the model, however, is restored if one in- 
cludes the quantum fluctuations; they are shown in Fig. 
|] as functions of the rotation angle 9. Note that the 
quantum corrections are small and comparable to those 
of the 3D HAF. They do depend on the rotation angle, 
as the orbital- wave dispersion does. If the dispersion is 
completely 2D-like, (0 = 7r/4) = 3J-^/l ±711 (k), cor- 
responding to 0^(0 = 0) j35|), the quantum corrections 
are largest, as this dispersion has a line of nodes along 
the r — Z direction, i.e., U)7 Q * = for < q < ir. 

It is interesting to note that this fully 2D orbital-wave 
dispersion leads to an orbital ordered state at zero tem- 
perature only. At any finite temperature fluctuations de- 
stroy the long-range order (as in the 2D HAF). In this 
way one finds an orbital-liquid state at non-zero temper- 
ature, similar to that obtained in a Schwinger-boson ap- 
proach tw Ishihara, Yamanaka and Nagaosa for the doped 
system.Ea In higher order spin-wave theory, however, it 
might very well be that a gap occurs in the excitation 
spectrum. We expect, however, that this gap, if it arises, 
is small, with its size being proportional to the size of 
quantum fluctuations. 

The cubic symmetry of the model at E z = is satis- 
fied by the quantum-corrected ground-state energy (Fig. 

, as the energy is invariant under rotations by any an- 
gle 9 — mr/6, where n is any natural number. This can 
be expected as such a rotation corresponds to merely 
directing the orbitals along different cubic axes and/or 
interchanging the sublattices A and B. 

At finite values of E z the spectrum becomes gapped 
(Fig. H) because the degeneracy of the classical ground 
state is removed. Note that this is different from a HAF 
in a finite magnetic field, because there one still has full 
continuous symmetry in the spin-flop phase with respect 
to rotations along the field axis, which justifies the exis- 
tence of a Goldstone mode. For small orbital fields the 
gap firstincreases up to a maximum reached at E z = 3 J 
(2 J), while for larger values the gap decreases and van- 
ishes at E z = 6J (4J) for a 3D (2D) model. This 
can be understood by making the analogy with an AF 
Ising model in a magnetic field. As the field increases, 
the energy of a spin-flip excitation decreases as the loss 
of interaction energy associated with this excitation is 
partly compensated by the energy gain due to the paral- 



lel alignment of the excited spin to the field. Approaching 
the phase transition from the spin-flop to FM ordering, 
this mode becomes softer and is found atzero energy ex- 
actly at the transition point. The increase of the gap for 
E z > 6 J is therefore due to a complete FM alignment of 
the orbitals along the field axis in this parameter regime, 
and any excitation acts then against the orbital field. In 
Fig. H also the quantum corrections to the ground-state 
energy and the order parameter are shown as functions of 
the orbital field. As expected, the presence of a gap in the 
excitation spectrum suppresses quantum fluctuations. 

The quasi-2D state, i.e., the alternation of \x) and \z) 
on two sublattices, is destabilized by any finite orbital 
field E z 7^ 0. Taking for instance E z > 0, the orbital 
field shifts downwards the energy of the occupied \x) or- 
bitals on one sublattice by E z /2, while the energy of \z) 
orbitals increases by E z /2 on the other sublattice, so that 
the orbital field does not lower the energy of the system. 
Therefore, the state with 9 = 7r/4 is selected instead fp*j 
the classical ground state for any non-zero orbital fieldE2l 
because in this state the pseudospins can tilt towards 
the orbital field, reducing the energy of the system. Such 
orbitals forming the ground state are shown in Fig. 0. 

The dispersion for the unrotated state suggests that 
the effective dimensionality of the system is reduced from 
three to effectively two dimensions in the absence of 
an orbital splitting. This anisotropy results in overall 
smaller quantum fluctuations in a 3D model at E z = 
than in the corresponding HAF, as discussed above. The 
anisotropy of the 3D model manifests itself in the energy 
contribution coming from the bonds in the (a, b) plane 
and along the c-axis, as shown in Fig. |t]. For E z = 
the classical energy along the chain is zero, and the or- 
bital superexchange energy of 3 J/4 is gained by the pla- 
nar bonds. One finds that quantum contributions to the 
ground-state energy tend to decrease the energy stored 
in the bonds along the chain, and increase the energy in 
the plane. 

For the 2D system the situation at orbital degeneracy 
is quite different (see Fig. ||). The lack of interactions 
along the c-axis breaks the symmetry of the model al- 
ready at E z = 0, opens a gap in the excitation spectrum 
and suppresses quantum fluctuations. Increasing the or- 
bital field, however, the system resembles the behavior 
of a 3D system, where also the gap closes at an orbital 
field which compensates the energy loss due to the or- 
bital superexchange between identical (FM) orbitals. At 
this value of the field (E z — 4 J), the full dispersion of 
the orbital waves is recovered, see Fig. ||, the spectrum is 
gapless, and the quantum fluctuations reach a maximal 
value. We would like to emphasize that this behavior is 
qualitatively different from the HAF, where the anoma- 
lous terms cx T^T^ and T^TJ are absent, and quantum 
fluctuations vanish at the crossover from the spin-flop to 
FM phase. 
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V. COMPARISON WITH EXACT 
DIAGONALIZATION IN A 2D CLUSTER 

It is instructive to make a comparison between the an- 
alytic approximations of Sec. Ill and the exact diagonal- 



ization of 2D finite clusters using the finite temperature 
diagonalization methodic In contrast to the mean-field 
approach presented in Sec. |n|, one finds a unique ground 
state at E z = by exact diagonalization, and no rotation 
of the basis has to be made to investigate the stability 
of the ground state. Making the rotation of the basis 
(fjj) is however still useful in exact diagonalization as it 
gives more physical insight into the obtained correlation 
functions which become simpler and more transparent 
when calculated within an optimized basis. Moreover, 
they offer a simple tool to compare the results obtained 
by exact diagonalization with those of the analytic ap- 
proach, presented in Sec. IV. We shall present below the 



results obtained with 4x4 clusters; similar results were 
also found for 10-site clusters. 

First we calculate nearest-neighbor correlation func- 
tion in the ground state (TfT* +R ), where the operators 
with a tilde refer to a rotated basis, 



T* = cos 2077 +sin2#Tf 
f z l+R = cos 2^T? +R + sin 2^Tf +R , 



(36) 



so that the correlation function depends on two angles: 
4> and ip. 

In Fig. |^ the intersite orbital correlation in the ground 
state is shown as a contour-plot. White (dark) re- 
gions correspond to positive (negative) values, respec- 
tively. One finds that the neighbor correlations have their 
largest value if the orbitals are rotated by cf> — 7r/4 and 
-0 = 37r/4 (or <\> = 37r/4 and ip = 7r/4), i.e., under this 
rotation of the basis states the system looks like a ferro- 
magnet, indicating that the occupied (\x) + |z))/V2 and 
(\x) — \z))/y/2 orbitals are alternating in a 2D model, as 
shown in Fig. @. Note that quantum fluctuations are 
small as in the ground state and one finds (T*T?, R ) ~ 
0.246. This demonstrates at the same time the advan- 
tage of basis rotation in the exact diagonalization study, 
because in the original unrotated basis one finds instead 
(TfT?, R ) ~ 0, which might lead in a naive interpretation 
to a large overestimation of quantum fluctuations. 

The ground-state correlations {T*T* +R ) are in excel- 
lent agreement with the LSWT results. By comparing 
the results obtained at T = 0.1J, 0.2 J and 0.5 J we found 
that the calculated low-temperature correlation functions 
are almost identical in this temperature range, and thus 
the values shown in Fig. |^ for T = 0.1J are represen- 
tative for the ground state. They demonstrate an insta- 
bility of the system towards the symmetry-broken state. 
This finding is similar to the HAF, where such a ten- 
dency could also be found, in spite of the relatively small 
size of the considered clusters £3 The intersite correla- 
tions decrease at higher temperatures, and thus the re- 
sults at higher temperatures T > 2J differ dramatically 



from those presented in Fig. ||. By investigating the tem- 
perature dependence of the orbital correlation functions 
it has been recently established that the orbital order 
melts at T ~ J.E3 

Next we discuss the results for the dynamical orbital 
response functions in the case of orbital degeneracy. The 
dynamic structure factors for the orbital excitations are 
defined as, 



7q Z » = «- / dt (2^T* q (t)) expt-iwt), (37) 

L ~ J— oo 

1 f°° 

7q + ~M = 7T / dt < T q +T -qW> exp(-iwt). (38) 



The corresponding orbital response functions evaluated 
with respect to the rotated local quantization axis, 
Tq Z (w) and (to), are defined instead by tilded op- 
erators ([36]). The LSW approximation does not allow 
to investigate the consequences of the coupling of single 
excitonic excitations to the order parameter, represented 
by the terms cx TfT?, as these terms contribute only in 
cubic order when the expansion to the bosonic operators 
is made. Therefore, the mixed terms involving products 
of T? and T? at two neighboring sites in the Hamiltonian 
Eq. (^) could contribute only in higher order spin-wave 
theory. This motivated us to perform similar calculations 
within finite temperature diagonalization for the simpli- 
fied Hamiltonian without these terms, which we refer to 
as the truncated Hamiltonian. 

The transverse and longitudinal orbital response func- 
tions, (u>) and T zz (ui), calculated for the full orbital 
Hamiltonian given in Eq. (j^), and with respect to the 
original choice of orbital basis {\z), \x)} are shown in 
Fig. || The response functions consist in each case of 
a dominant pole accompanied by a pronounced satellite 
structure, both for the transverse and longitudinal exci- 
tations. The dominant pole energies are close to the two 
modes found by the spin-wave analysis, but the energies 
are lowered for those momenta which show strong satel- 
lite structures at higher energies. In order to establish 
that these satellite structures are, in spin-wave language, 
due to higher-order processes in spin- wave theory, we per- 
formed the same calculation for the truncated Hamilto- 
nian. One finds then an excellent agreement (Fig. |^) 
between the single modes which result from the numer- 
ical diagonalization and the position of the orbital-wave 
as determined by Eq. ([33]). 

An additional numerical check was made using the ro- 
tated orbital basis which determines the transverse cor- 
relations (T+TTq) shown in Fig. Here (j> — ir/4 
and ifj = 37T/4 for A and B sublattice, respectively. In 
contrast to the unrotated orbital basis, we find that the 
dispersive poles in the longitudinal excitation spectrum 
vanish entirely. Instead, this missing mode appears now 
in the transverse channel, where both modes are found 
and are accompanied by the satellite structures. We note 



7 



that the response functions in the rotated basis resem- 
ble the response functions of a HAF, where the spin- 
wave excitations ace also found only in the transverse 
response function.E3 As before, the satellite structures 
disappear and two distinct peaks with the same inten- 
sities are found, if the truncated Hamiltonian is used. 
These two peaks merge into a single structure at those 
values of q at which uj^(<f) = 0) = tJ^((j> = 0). 

The successful comparison of the LSWT results with 
the numerical diagonalization is summarized in Fig. [ll], 
where the first moment of the spectra of the full Hamil- 
tonian, the single mode which results from the trun- 
cated Hamiltonian, and the dispersion from the LSWT 
[as found from Eq. (|3^) at (j> — and neglecting the 
dispersion due to fcJ, 



± = 3jJl±± 7 ||(q), 



(39) 



are compared. Clearly, Fig. [ll] shows an excellent agree- 
ment between the dispersions calculated within different 
approaches. From the fact that the first moment of the 
response functions for the full Hamiltonian coincides with 
the response functions of the truncated Hamiltonian, we 
conclude the LSWT approach captures the leading terms 
in orbital dynamics, and the satellite structure for the full 
Hamiltonian is due to higher-order orbital-wave interac- 
tions that are not taken into account within LSWT. 

As already discussed in Sees. || and [V, the orbitals 
are tilded in the direction of the effective orbital field due 
to the orbital splitting. This behavior is also found in our 
diagonalization studies, and one finds that the expecta- 
tion values (TfT?) and (TfT/) depend on the orbital 
splitting E z (Fig. |l2|). Both correlation functions are 
calculated using the unrotated basis-set. Although the 
crossover to the uniform state with |a;) orbitals occupied 
is smooth in the numerical study, the agreement with the 
result of the LWST approach is once again very good. 



VI. SUMMARY AND CONCLUSIONS 

In summary we considered the low-energy effective 
model for strongly correlated electrons in two-fold degen- 
erate e g bands and investigated the consequences of the 
orbital degrees of freedom for a FM system at half-filling. 
Although a saturated ferromagnet might only be realized 
in undoped systems in high magnetic fields, we believe 
that this model is a relevant starting point in order to un- 
derstand the influence of orbital dynamics in FM doped 
transition metal oxides, such as the colossal magnetore- 
sistance manganites. The FM state can be considered as 
the reference state into which holes are doped. 

We have shown that the orbitals order and alternate 
between the two sublattices, which is equivalent to the 
AF order in the (orbital) pseudospin space. The or- 
bital excitation spectrum consists of two branches, one 
of which is found to be gapless in the 3D system within 



LSWT. This is a consequence of the degeneracy of the 
model at a classical level. Thus, we find here a quite 
peculiar situation. Although the Hamiltonian is not in- 
variant under a global orbital rotation 9 [see Eq. (7)], 
and the correlation functions along the different direc- 
tions change with 9, this rotation does not affect the 
energy of the classical ground state. It turns out, how- 
ever, that quantum fluctuation corrections to the ground 
state energy restore the cubic symmetry and select the 
rotation angle at either 9 = 0, or 7r/6, or 7r/3. For these 
states the orbital excitation spectra are purely 2D, which 
demonstrates that in fact the effective dimensionality of 
the orbital model ^) is reduced to two. We note that for 
a gapless 2D dispersion of the orbital excitations in three 
dimensions one expects that orbital long-range order dis- 
appears away from T = 0, like in the 2D Heisenberg 
model. 

The orbital model (|l|) has novel and interesting quan- 
tum properties, and is quite different from the Heisenberg 
model. In contrast to the HAF, the quantum fluctuations 
in zero field are smaller in the 2D than in the 3D orbital 
model, and the 2D model is therefore more classical due 
to the anisotropy of the orbital interactions. However, a 
finite splitting of the orbitals E z opens a gap in the exci- 
tation spectrum of a 3D model, and the quantum effects 
in an orbital-flop phase stable at small values of \E Z \ are 
then reduced. Both in a 3D and 2D case, the orbital ex- 
citation spectra become gapless at the crossover value of 
E z at which the 'FM' orbital order sets in. The quantum 
fluctuations are then identical and reach their maximum. 

We have verified that the linear approximation within 
the spin- wave theory reproduces the essential features of 
the orbital excitation spectra in the present situation, 
where the total pseudospin quantum number and the zth 
component of total pseudospin are not conserved. First 
of all, the results of the LSWT and the first moments 
of the orbital excitation spectra obtained from the finite 
temperature diagonalization yield two modes in excellent 
quantitative agreement in a 2D model, when the anoma- 
lous interactions ~ TfT* are neglected. Second, the in- 
clusion of these terms in exact diagonalization leads to 
a satellite structure in the response functions, while the 
first moment remains almost unchanged. In order to cal- 
culate the satellites observed in the exact diagonalization 
one needs to go beyond the leading order in the spin- 
wave theory, and include the cubic terms in the boson 
operators. The comparison with finite temperature diag- 
onalization allows us to conclude that orbital waves with 
the dispersion of the order of 4 J are characteristicjor the 
FM states of the undoped degenerate e g systems^ 

We believe that the presented analytic treatment of 
orbital excitations provides a good starting point for 
considering the dynamics of a single hole in the orbital 
model. This would allow to understand the origin of 
puzzling optical and transport properties of the doped 
manganites and clarify the expected differences with the 
hole dynamics in the t — J model. 
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TABLE I. Occupied orbitals at two sublattices A and B as found in the classical FM ground state at orbital degeneracy 
(E z — 0) for a few representative values of the orbital rotation angle 9 [see Eq. (ffl)]. 
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tion of the virtual d t dj —* cqdj excitations in LaMnC>3 for 
the starting FM configuration d\ z ^d^ x ^\Q) which involve the 
high-spin \ a A\) state and generate effective orbital superex- 
change interactions: (a) for a bond along the c-axis, (ij) _L; 
(b) for a bond within the (a, 6) plane, (ij) \\. 




FIG. 2. Alternating orbital order in FM LaMn0 3 : 
{\x) + \z))/V2 an - \z))/V2 as found at E z -> 0. 
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FIG. 3. Orbital-wave excitations as obtained for different 
values of the crystal-field splitting E z for a 3D (left) and 2D 
(right) orbital superexchange model (^). The result shown for 
a 3D system at E z — was obtained for the orbitals rotated 
by 9 = 7r/4 (Q), and corresponds to the E z — > limit of the 
orbital-flop phase. 



0.00 



< 
A 

I 

I 



-0.02 



-0.04 



-0.06 



-0.08 



-0.10 











AT Z 






AE 












/ 


/""\ 




/ \ 










Heisenberg 3D -- 



7[/6 71/3 

Rotation angle 6 



71/2 



FIG. 4. Quantum corrections for the 3D system as func- 
tions of rotation angle 9 for: the renormalized order parameter 
AT Z (full lines), and the ground-state energy AE/J (dashed 
lines). 
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FIG. 8. Contourplot of the rotated nearest-neighbor or- 
bital correlation function {T Z T Z +R ) as function of the angles <j> 
and ip for a 16-site planar cluster with E z — and T = 0.1J. 
White regions correspond to positive (FM) and black areas 
to negative (AF) orbital correlations, i.e., (T Z T Z +R ) > 0.24 
(< —0.24), respectively. They are separated by 25 contour 
lines chosen with the step of 0.02 in the intervall [-0.24,0.24]. 



3 



A 

v 

— 

20 





Ji/V^ 


- q=(3i,ic) ft 




- q=(jc,n/2) » 

J Si 




- q-M) 

A I 


V_ 


- q={Ji/2,0) 

A 




q.(0,0) 


H M1 - 



q=(7C/2,7C/2) 




J 




q={7c,jt/2) 




q=(jC,0) 

a! 




q=(7t/2.0) 

A i' 




q=(0,0) 

A' 


H full 
U 

"trunc 



FIG. 9. Transverse response function Tq (w) (left part) 
and longitudinal response function T^ z (uj) (right part) calcu- 
lated with respect to the original basis states (9 — 0) using 
finite temperature diagonalization at T = 0.1J for a 16-site 
planar cluster with E z — for: (a) the orbital interactions 
as given by Eq. ( p"o"| ) (full lines); (b) the truncated Hamil- 
tonian (dashed lines). Vertical lines indicate the position of 
orbital-wave excitations as obtained from LSWT. The spectra 
are broadened by V = 0.1J. 
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FIG. 11. Dispersion of orbital waves along the main direc- 
tions in the 2D Brillouin zone calculated for a 4 x 4 cluster. 
Results for the first moment calculated using the full and 
truncated Hamiltonian (empty and full symbols) are com- 
pared with the dispersions for the two modes as obtained by 
the LSWT, shown by a solid and dashed line, respectively. 



FIG. 10. Transverse response function Tq (w) for the ro- 
tated orbitals as a function of frequency uj for different mo- 
menta at low temperature. Calculations were performed for 
a 16-site 2D cluster at E z = and T = 0.1J for: (a) the 
orbital interaction as specified in Eq. ( Jlo| ) (full lines), and 
(b) neglecting the mixed terms oc T^Tj . The spectra are 
broadened by T = 0.1J. 




FIG. 12. Ground state expectation values (TfT/) (left) 
and (TfT/) (right) for a 16-site 2D cluster (E z = 0, T = 0.1J) 
as functions of the orbital splitting E z /J (squares). Full 
lines represent the same expectation values obtained using 
the LSWT approach. 
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